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We study Skyrmion quantum numbers, charge and statistics, in (2 + 1) dimension induced by 
quadratic band toucing(QBT) fermions. It is shown that induced charge of Skyrmions is twice 
bigger than corresponding Dirac particles' and their statistics are always bosonic. Applying to the 
Bernal stacking bi-layer graphene, we show that Skyrmions of quantum spin Hall(QSH) are charge 4e 
bosons, so their condensation realizes charge 4e superconductivity(SC). The phase transition could 
be a second order, and one candidate theory of the transition is 0(5) non linear sigma model(NLSM) 
with non-zero Wess-Zumino-Witten(WZW) term. We calculate renormalization group beta function 
of the model perturbatively and propose a possible phase diagram. We also discuss how QBT 
fermions are different from two copies of Dirac particles. 



Introduction : Electronic strucutres with helicity such 
as ones of Dirac particles and quadratic band touch- 
ing(QBT) fermions in two-spatial dimension receive a lot 
of attention especially with the discovery of graphene^ 
The helicity and its associated Berry phase play an im- 
portant role in quantum mechancal phenome na such as 
half-integer Quantum Hall effects in graphene P^l In par- 
ticular, QBT fermions lead more interesting and puz- 
zling physics due to a stronger interaction effect and 
are not well understood yetP"^ Not only graphene but 
also some spin liquids and topological insulator also have 
QBT fermions structures EES! 

In (2+1) dimension, there is another common and 
interesting excitation, Skyrmion of 0(3) vector order 
parameters ! 11 * 12 ! For example, the Neel order phase 
has Skyrmions because its ground state manifold is 
a two sphere due to its broken symmetry structure, 
SU(2)/U(1) — S 2 . Basically, a Skyrmion number is a 
wrapping number of a given configuration on a two di- 
mensional space with identified boundaries. 

One natural question with the two excitations, QBT 
fermion and Skyrmion, is how they interact with each 
other. In this paper we answer the question focusing on 
Skyrmions' induced quantum numbers by QBT fermions. 
Extensiv e rese arch already exists on Skymions with Dirac 
particlesPEsl 

Our goal is two-fold. First, we show Skyrmions receive 
quantum numbers by interacting with QBT fermions, 
which are doubled compared to Dirac particles'. In ad- 
dition, we apply the mechanism of the induced quan- 
tum numbers of Skyrmions to the Bernal stacking bi- 
layer graphene structure and study an exotic quantum 
phase transition in the bi-layer graphene. We report two 
findings. First, quantum spin hall phase in the bi-layer 
graphene is connected to charge 4e SC by QSH Skyrmions 
condensation. And second, the phase transiiton could 
be a second order and we propose one candidate the- 
ory, 0(5) non linear sigma model with the Wess-Zumino- 
Witten term. 

In terms of charge 4e SC, one can see that it is almost 
impossible to realize it with the conventional BCS theory 
because two pairs of charge 2e order parameters have a 



better chance of forming than charge 4e SC. Therefore, 
one key factor to achieve charge 4e SC is how to prohibit 
charge 2e order parameters. One possibility of charge 4e 
superconductivity was suggested by Berg, et, aZ.P^Main 
idea of their study is to prohibit 2e pairings thermally 
to make 4e superconductors, so they only find non-zero 
temperature charge 4e SC. Our research is a significant 
departure from the previous one by Berg in that our 4e 
SC is a ground state at zero temperature. 

Abanov and Wiegmann studied a non-BCS way to 
achieve superconductivity focusing on acquired quantum 
numbers of topological objects interacting with Dirac 
fermionsPS If a quantum number-endowed topological 
object is a boson, it is possible to condense the object, 
and it induces quantum phase transition as a textbook 
example in the two dimensional XY model. Applying the 
idea to monolayer graphene, Grover and Senthil argued 
that quantum spin hall phase in monolayer graphene can 
be connected to charge 2e superconductor through the 
CP{\) deconfined quantum criticalityPil Below we show 
how QBT fermions change Skyrmions' properties and as- 
socited Quantum phase transition. 

Model Hamiltonian : Before going into QBT fe rmoins, 
let us briefly review results of Dirac particles ! 14 * 15 * 21 ! 
Hamiltonian consists of Dirac fermions with 0(3) vector 
order parameter. 

H D irac = ^ [Ko x + k y a v + ga z n ■ f] ip. (1) 

, where a 1 , t j are different kinds of Pauli matrix, so it is 
necessary to introduce at least four component spinors. It 
is well known that the 0(3) vector fields allow non-trivial 
topological configuration, Skyrmion, due to ^(S 2 ) = Z. 
Its topological current is = ^e^xri ■ [d v n x d\n) 
and charge is spatial integration of its time component, 
gSkyr = j d 2 x±n ■ (d x n x d y n). 

It was shown that the induced fermionic current, < 
jf erm > ; i s proportional to Skyrmion current and its 
charge is exactly determined by Skyrmions number upto 
number of flavors. 

Q/erm = Nf QSkyr ^ 



Also, it was well known that Skyrmions' statitical prop- 
erties could be read off by examining two Skyrmions 
exchange property or measureing a phase of 2tt rota- 
tion. Those information is contained in the so-called Hopf 
term, which corresponds to ^(S* 2 ) = Z mathematically 
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, where Nf is the number of flavors and the gauge 
field(v4 M ) is defined by the relation, J M = s^ u \dvA\. 

Massage of the above discussion is manifest. 
Skyrmions of order paprameters receives charge and 
statistics by Dirac fermions. Dirac particles with even 
number flavors make Skyrmion bosonic charged parti- 
cles and odd number's make it fermionic charged par- 
icles. Notice that there is one dimensionful coupling 
constant (g) but it does not appear in our final results 
except its sign. It is a signal of topological origin of in- 
duced quantum numbers of Skyrmions. 

Now, let us consider quadratic band touching fermions 
interacting with 0(3) vector. In general, Hamiltonian is 



H = ^ 



(k\ - k 2 y )a x + (2k x k y )gy 
2m 



+ ga n ■ f 



This is a generalization of the Eqn. [T] from the p wave 
helicity to the d wave helicity with introducing one more 
dimensionful scale, effective mass(m). Obviously, the 
Lorentz symmetry is broken by quadratic dispersion re- 
lation. Such systems can be realized, for example, in 
bi-layer graphene or in the hole-doped semiconductors 
with the inversion symmetry. 

One way to see induced charge of Skyrmion is to eval- 
uate fermion's energy spectrums under Skyrimon back- 
ground. Turning on Skyrmion configuration adiabati- 
cally, one can show there are zero energy level crossing 
states. For one spatial dimension with 0{2) symmetry, 
Yao and Lee showed there were two level crossing for the 
QBTP In principle, one can do similar calculation in 
two spatial dimension. Instead, we use a field theoritic 
method, the graident expansion, to see both Skyrmion's 
charge and statistics also. Similar calculation method has 
been applied to Dirac particles in various literatures.^ 

For QBT fermions, broken Lorentz symmetry makes 
the expansion little more tedious but it is straightforward 
to do if one only keeps the lowest term. Partition function 
is 



Tr(e" 



- J d 2 xdrC 



, where the Lagrangian is 

C = ^Dxjj = ip^d T ip - 



H. 



(3) 



One key step of the gradient expansion is to consider 
variation of the effective action. 
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FIG. 1. Induced current(a) and Hopf term(b) diagram. The 
plain line is for the propagator. (Go) The square vertex is for 
the variation of the D operator. (SDD^) The circle vertex is 
the coupling between the fermion and the order parameter 
variations. (M) 



Following commutation relations are useful. 

[kl - k 2 y ,h] = 2{-i){k X d X h - kyOyh) + {dl - d 2 y )h 

[k x k y + k y k x , n\ — 2(—i)(k x d y h + k y d x n) + 2d x d y h 
If we keep only the lowest order terms, we get 

k 2 
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With these expansion, one can do perturbative calcula- 
tion with Green's function, Go with perturbation, M. In 
Fig. [Tj we illustrate two diagrams for induced current (a) 
and Hopf term(b) calculations. Note that to calculate 
induced currents it is necessary to couple external gauge 



field as in an usual way, {k^ 
the first diagram, we find 



— > fc M — a M ). By evaluating 



2N f Q skr . (5) 

Here, we put Nf back to the formular for notational pur- 
pose. 

Few remarks are in order. First, there is a factor two 
compared with the Dirac fermion's of Eqn. [2j The in- 
duced charge of the Skyrmion by QBT fermions is dou- 
bled compared to the Dirac case. Second, the result in- 
dicates that Lorentz symmetry is not a crucial factor for 
topologically induced currents. In other words, the vac- 
uum of the non relativistic fermion with the topologi- 
cal background could be polarized and its polarization is 
twice bigger. Also, one notice that there is no effective 
mass dependence in the current calculations. This fact 
is another signature of its topological origin. Otherwise, 
one expect that the effective mass term would appear in 
momentum integration. 

To complete induced quantum numbers, let us study 
statistics of topological defects. It was shown by Wilczek 
and Zee that the Hopf term of the non-linear sigma 
model, so-called the "theta" term, describes the statistics 
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of the Skyrmion. 
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, where 9 term is determined by microscopic models. If 
we consider the rotation(or exchange) of the soliton, then 
its current naturally induces the non-trivial imaginary 
factor depending on its spin(J). 



S^S + ie (S^S + i2irJ) 



(7) 



So the Skyrmion spin is J = and tt(2tt) corresponds 
to fermion(boson). Other values indicate anyonic statis- 
tics of the Skyrmion. 

Following the same steps as in the current calculaiton, 
it is straighforward to obatin $S q °*J using the second 
diagram of Fig. [ljb). We find that imaginary part of ef- 
fective action's variation corresponding to the Hopf term 
is 
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, where h = f • ft is used. 

It is obvious that the variation becomes zero, which 
naively looks no Hopf-term contribution. However, it 
is well known that the Hopf term caluculaiton is more 
subtle due to its non-trivail homotopy group, ^(S 2 ) = 
Z. In other words, the zero variation is an artifact of our 
perturbative method. To ove rcome this difficulty, we use 
the embedding method? 15 * 23 ! Main idea is to extend the 
order parameter group, CP 1 
avoid the ambiguity. 



S 2 , to enlarged CP M to 



n ■ t — > 



2zz f - 1 



(9) 



with z t — (zi, • • • , zm+i) and z^ z = 1. Due to the fact, 
n^(CP M ) = 0, for M > 1 it is safe to use a perturbation 
theory. Variation of the effective action has no ambiguity 
and our perturbative calculation is well-defined. Then, 
we take the limit M — > 1 and read off the induced Hopf 
term. And the induced Hopf term is 



(10) 



We put Nf back to the final forms for notational pur- 
poses. As we can see, the factor two again appears in 
the Hopf term calculation, which is consistent with the 
charge current calculation. Therefore, Skyrmions with 
the QBT fermions are always bosons and their charge 
is determined by number of flavors(2iV/). Again, the 
dimensionful effective mass does not appear due to its 
topological origin. 

How do we understand Skyrmion charge doubling ef- 
fect in QBT fermions? One way to understand is to add 
small perturabation which do not break other symme- 
tries such as gauge potential disorder with strengh(q). 
Then, it is easy to see the QBT point at the minium 
splits into two Dirac points below the energy sclae cor- 



responding (q). If only low energy properties are dom- 
inant to determine topological charges, we can just add 
two Dirac particles to make QBT fermions, which ex- 
plains the doubling effect. It seems the argument that 
QBT is equivalent to two Dirac particles is compelling, 
but they show different fermionic mass structures. Espe- 
cially, competinig order parameter descriptions based on 
the mass term structure are different! 3 -^ 

Physical Realization : One of the well-known examples 
of the QBT fermion is a bi-layer graphene with Bernal 
stacking. Also, hole-doped GaAs with the inversion sym- 
metry along the perpendicular axis has similar structure. 
In this paper, let us focus on the bi-layer graphene. Its 
tight binding Hamiltonian is 



H = -t i c \,i c hi + 4,j c 2,i) + * ( c m c 2,* + h.c.) 

(11) 



<i,j> 



iGA-B 



, where c\^) % ^ s a creation operator of the first (second) 
layer at the site(i), and the second summation is only for 
one sublattice where sites on one layer is on top of sites 
on the other layer, {A — B). 

Among four band dispersions, only two bands become 
important and the low energy Hamiltonian becomes 



{2k x k y )o-y p z 



2m 



V; (12) 



, where a, t, p correspond to sub-lattice(or layer), spin, 
and valley Pauli matrix. This Hamiltonian has the helical 
quadratic dispersion structure as our model Hamiltonian. 
Only considering non-valley mixing states, there are two 
order parameters with 0(3) vector properties. 



Oqsh = V'Vfp 2 ?/' 
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The first one corresponds to the Quantum Spin 
Hall(QSH) order parameter because it does not break 
the time reversal symmetry even though the spin rota- 
tional symmetry is broken. The second one breaks both 
the spin rotational symmetry and time-reversal symme- 
try, so it's a Neel order parameter. To make connnection 
to our model Hamilnotian manifest, let us rotate the ba- 
sis as follows. 



i 



1 - a x 



P z )i> 



(14) 



Then the p dependence of the non-interacting Hamilto- 
nian i?G,o dropped. Of course, the order parameter terms 
also rotate under the transformation and become 
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l/j^ (J Z Tip , N ->^<T Z T P Z ll> (15) 



Therefore, the QSH order parameter exactly corresponds 
to our model Hamiltonain except doubling of the number 
of flavors (valley index). 

Now we can apply our model Hamiltonian results di- 
rectly. Skyrmion of the QSH order parameter is a boson 
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with charge four considering the valley index. On the 
other hand, if we consider the Neel order parameter, it 
is easy to see the fermionic charge is zero and the Hopf 
term becomes trivial because the signs of the coupling 
constant of different valley fermions are opposite and 
their contributions are cancelled out. It indicates that 
the Skyrmion of the Neel order parameter is a charge 
neutral boson. Therefore, it is manifest that Skyrmion's 
induced quantum number depends on order parameters' 
property, especially how an order parameter couples to 
Dirac particles with different hclicity. 

Let us focus on the QSH phase. By tuning interac- 
tions, one can destroy the QSH order parameter. Similar 
to vortex proliferation in XY model in (1+1) dimension, 
it is possible to eliminate the QSH phase by condensing 
the Skyrmion due to its bosonic properties. This phase 
cannot be the same as the semimetal phase because the 
Skyrmion's electric charge is four. Therefore one can re- 
alize the charge 4e superconductors by condensing topo- 
logical defects of the QSH order parameter. This scenario 
does not require any Fermi surfaces and the Cooper pair 
formation, so there is no way to form charge 2e order 
parameters here. In other words, we can circumvent the 
instability problem of the charge 4e SC into formation of 
two charge 2e Cooper pairs. Therefore, it is one concrete 
way to achieve charge 4e SC through the quantum phase 
transition from QSH. 

What kind of a quantum phase transition is the tran- 
sition from QSH to charge 4e SC? Because two phases 
break different symmetries, it should be a first order 
phase transition under the Landau-Ginzburg paradigm. 
But, a second order phase transition is also possible via a 
deconfined quantum criticality! 20 * 21 ! In quantum mechan- 
ics, topological defects of one order parameter can carry 
the other order parameter's quantum numberPSl Their 
condensation suppresses the original ordered phase and 
induces the different symmetry-broken phase simultane- 
ously. Thus, a continuous phase transition can naturally 
appears in this process. 

One candidate critical theory of the phase tran- 
sition is 0(5) model with the Wess-Zumino-Wittcn 
term(WZW) ! 24 l 27 l Introducing a multiplet, <f> — 
(ft, Re ip ie , Im^ 4e ), the model is 

S = [ drd 2 x ^(d$) 2 + 2mkT wzw ' ($) (16) 

J %g 

T wzw ($) = 

, where k is an integer number level. 

Without the WZW term, the beta function of g is well 
understood. There are two fixed points; g = corre- 
sponds to the ordered phase and g = g c describes a criti- 
cal fixed point. Beyond the critical fixed point, the beta 
function monotonically increases. One can expect the 
WZW term plays non-trivial effect due to its topological 
orgin. Indeed Witten showed the WZW term in (1 + 1) 



dimension induces stable fixed points beyond an unstable 
fixed point PU 

Here, we follow his method and do the beta function 
calculation perturbatively. We find that 

^=/3 (5 2 )"O fcV° + O(g 12 ,fc 4 ). (17) 

The /3q function is one without the WZW term, which is 
obtained from the literature!^ 

/3o(5 2 ) = -.9 2 +3.g 4 + 3g 6 + |5 8 + 365 10 

, where b = ^ + 3£(3) and dimensionless coupling con- 
stant g is introduced to match higher loop calculation. 

We note that our calculation is not controlled basi- 
cally because of existence of the original mass scale of 
the coupling constant. In contrast, 0(4) NLSM on 1 + 1 
dimension calculation is controlled by 1/k expansion, and 
one can show that there are stable fixed points with non- 
zero k values, the so-called SU(2) k WZW fixed point.^ 
One physic al realization of such fixed points is spin 1/2 
chain! 28 * 29 ! Even though our calculation is not controlled, 
it clearly shows that the topological term induces a non- 
trivial effect. The topological term's contribution is of 
the opposite sign to the other terms in the beta func- 
tion. Physically, this means that the topological term 
suppresses the tendency of becoming a disordered state. 

There are three possible scenarios based on our per- 
turbative calculation. First, the topological term con- 
tribution is not big enough, so it does not induce any 
stable fixed point. With a small enough value of k, the 
beta function is barely affected by the topological term. 
Phase transition is described by the 0(5) NLSM model 
without the WZW term and the topological contribution 
only gives a small correction. The second scenario is that 
the topological term is so big that the beta function be- 
comes always negative. An implication of this case is 
that the system is always in the ordered phase. Both 
cases can be described by the 0(5) NLSM without the 
WZW term, so there is no new physics. The final case 
is that the topological term and the higher order fluc- 
tuation of the order parameter have similar amounts of 
effects. Thus, one new stable fixed point appears beyond 
the usual critical fixed point. 

This new fixed point is our candidate for the deconfined 
quantum criticality between QSH and charge 4e SC. By 
introducing symmetry breaking term such as y(ft ■ ft — 
IV'ScI 2 ), each phase can be accessed from the stable fixed 
point. Numerically, we find Co = J+ and k c = 2 with a 
hard-cutoff regularization schemePS If one believes that 
the SU (2)k WZW fixed points in 1 + 1 dimension have its 
extension to ones in the 2 + 1 dimension, our new fixed 
point would be a natural candidate. Figure [2] illustrates 
both the three possible cases in terms of renormalization 
group flow and the phase diagram at the critical level. 

Finally, QBT fermions are different from two copies of 
Dirac fermions in terms of fermion mass term structures, 
which usually explain competing order parameters. As 
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FIG. 2. (Color Online) (a) Beta function behaviors of 0(5) 
NLSM model with different levels WZW level, k. The dot- 
ted(black), normal(red), and dashed(blue) lines correspond to 
fc c cases. At k — k c , there is one stable 
interacting fixed point(balck dot), (b) Phase diagram with 
RG flows for the critical level at k = k c . 

shown by HerbuP^ recently, there is a "hidden" internal 
structure(pseudo-spin) of competing order parameters of 
Dirac particles. For mono-layer graphene with Dirac par- 
ticles, SC order parameter can make competing 0(5) vec- 
tor with either Quantum spin Hall or density wave order 



parameters. In other words, there are two possible ways 
to make competing 0(5) order parameters. However, one 
can easily check that QBT fermions do not allow the same 
two possibilities. Only one choice is allowed. For the SC 
order parameter, only the QSH order parameter is pos- 
sible to make 0(5) vector. In that sense, QBT fermions 
describe a more restricted competing orders' structure, 
which originates from the difference of Clifford algebra 
representations of Dirac particles and QBT fermions.^ 
We leave this interesting problem for the future work. 

Conclusion: In this paper, we showed how Skyrmions 
interacting with QBT fermions obtains their quantum 
numbers. It is shown that QBT fermions double induced 
charges compared to Dirac fermions' and make all the 
Skyrmions bosonic. Applying the mechanism to bi-layer 
graphene, we also showed that QSH Skyrmions in bi- 
layer graphene are charge 4e bosons. By condensing the 
Skyrmions, it is possible to achieve charge 4e supercon- 
ductivity. Such mechanism is far different from the gen- 
eral BCS-like description, and we argue that it might 
be described by the deconfined quantum criticality. We 
propose the 0(5) non-linear sigma model with the WZW 
term at level k c as a critical theory. 
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